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ABSTRACT 

ARL  is  pursuing  the  goal  of  developing  a  finite- 
element-based  design  methodology  for  thoracic  body 
armor.  We  describe  progress  in  modeling  two  essential 
ingredients,  a  Kevlar  vest  and  the  human  thorax. 

1.  INTRODUCTION 

Thoracic  armor  is  generally  designed  by  an 
experimental  methodology,  which,  since  the  1980s,  has 
been  based  on  the  clay-backed  test  codified  in  an  NIJ 
standard  (National  Institute  of  Justice,  1987).  This  test 
involves  placing  the  armor  on  a  standardized  block  of 
clay,  shooting  the  armor  with  the  “design  threat,”  and 
then  measuring  the  depth  of  the  crater  left  in  the  clay.  No 
account  is  taken  of  hit  location  on  the  thorax,  and  as  a 
result  Kevlar  vests  are  currently  not  tapered  to  save  on 
weight.  Nor  is  account  taken  of  the  time  sequence  of 
thoracic  events  as  a  response  to  the  impact;  only  the 
permanent  deformation  in  the  clay  is  measured. 

Since  at  least  the  early  1990s,  the  automotive 
industry  has  been  engaged  in  the  development  of  3D 
finite  element  (FE)  models  of  the  human  thorax  for  use  in 
numerical  crash  simulations.  In  1999  ARL  explored  the 
potential  for  transitioning  to  the  field  of  body  armor 
design  these  advances  in  FE  thorax  modeling.  The 
Wayne  State  thorax  model  (WSTM)  (Wang,  1995)  was 
found  to  be  the  most  anatomically  detailed  model 
available.  Figure  1  shows  that  model’s  representation  of 
the  skeletal,  circulatory,  and  respiratory  systems. 

ARL  purchased  access  to  the  WSTM  in  2000  to  use  it 
as  a  starting  point  for  the  development  of  a  new, 
simulation-based  methodology  for  body  armor  design. 
Figure  2  shows  an  overview  of  our  proposed  approach. 
An  FE  model  of  a  body  armor  design  is  placed  on  the 
thorax  model,  and  the  two  are  then  impacted  with  a  model 
of  the  threat.  The  three  models  each  consist  of  an  FE 
mesh  (discretized  geometry/anatomy)  and  a 
representation  of  all  materials’  mechanical  properties. 
The  simulation  is  run  on  an  FE  code  such  as  LS-DYNA 
(Livermore  Software  Technology  Corp.,  2003).  The 
immediate  output  is  local  stresses  and  displacements 
throughout  the  thorax,  which  must  then  be  related  to 
quantitative  injury  assessment.  The  designer  then 
systematically  modifies  the  armor  and  observes  the 


resulting  changes  in  the  simulation  results,  thereby 
seeking  an  optimal  design. 


Fig.l.  Skeletal,  cardiovascular,  and  respiratory  systems 
as  represented  in  the  Wayne  State  Thorax  Model. 


_ _ 
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Fig.  2.  Calculation  overview. 

In  Section  2  the  WSTM  is  applied  to  the  case  of  the 
M882  bullet  at  445  m/s  versus  a  multi-ply  Kevlar  vest 
plus  thorax.  For  this  situation,  accelerometer  data  from 
instrumented  human  thoracic  tissue  were  available  to 
serve  as  benchmarks  for  the  simulations  (Mackiewicz  et 
al.,  in  preparation). 

Section  3  describes  an  orthotropic,  nonlinearly  elastic 
constitutive  model  for  a  Kevlar  vest.  Section  4  describes 
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progress  towards  the  development  of  a  thorax  FE  model 
that  is  computationally  more  robust  than  the  WSTM. 


2.  AN  APPLICATION  OF  WSTM  TO  A 
BALLISTICS  SIMULATION 

NIJ  Level  3A  body  armor  is  intended  to  offer 
protection  against  semiautomatic  weapons  and  generally 
takes  the  form  of  a  fabric  vest.  National  Institute  of 
Justice  (1987)  specifies  the  Level  3  A  “design  threat”  to  be 
the  M882  bullet  impacting  at  a  minimum  velocity  of 
427  m/s.  This  bullet  is  a  standard  NATO  round 
composed  of  a  lead  core  and  a  copper  gilding  jacket.  Its 
diameter  is  9  mm,  and  its  mass  is  8  grams  (Fig.  3). 


In  2000  tests  were  performed  at  the  Armed  Forces 
Institute  of  Pathology  (AFIP),  in  which  an  M882  was 
fired  at  445  m/s  into  a  multi-ply  Kevlar  vest  worn  by 
human  thoracic  tissue.  The  vest  was  composed  of  plain- 
woven,  600-denier  Kevlar  KM2.  The  bullet  was  aimed  at 
the  center  of  the  sternum.  Endevco  accelerometers  were 
surgically  implanted  at  four  locations  in  the  thoracic 
tissue:  the  posterior  sternum  (Fig.  4a),  the  spinuous 
process  of  the  T7  vertebra  (Fig.  4b),  the  carina 
(bifurcation  point)  of  the  trachea,  and  the  ligamentum 
arteriosum  (Mackiewicz,  in  preparation). 

The  AFIP  tests  were  simulated  using  the  LS-DYNA 
FE  code  (Livermore  Software  Technology  Corp.,  2003.) 
An  FE  model  of  the  vest  was  added  to  the  exterior  of 
WSTM  using  HyperMesh  software  (Altair  Engineering, 
2000)  (Fig.  5).  The  vest’s  mesh  consisted  of  one  8-node 
solid  (“brick”)  element  through  the  thickness.  The  typical 
vest  element  had  in-plane  edge  lengths  of  about  10  mm 
(consistent  with  the  typical  element  size  of  WSTM)  and  a 
3.63-mm  thickness.  This  thickness  is  the  ratio  of  the 
vest’s  known  areal  density  and  the  volumetric  density  of 
1440  kg/m^  that  has  been  used  for  Kevlar  (Johnson,  et  al., 
1999).  Isotropic  Hooke’s  law  was  applied,  with  a 
Young’s  modulus  of  74  GPa  and  a  Poisson’s  ratio  of  0.2. 


These  values  were  also  obtained  from  Johnson  et  al. 
(1999).  The  M882  bullet  was  modeled  with  a  single  lead 
element,  to  which  perfect  plasticity  was  applied  with  an 
initial  density  of  1 1350  kg/m^,  an  elastic  shear  modulus  of 
5.52  GPa,  and  a  flow  stress  of  34.5  MPa. 


Fig.4.  Two  accelerometer  locations  in  the  AFIP  tests:  (a) 
posterior  sternum,  (b)  spinuous  process  of  T7  vertebra. 


Fig.  5.  Kevlar  vest  added  to  the  WSTM. 

Figure  6  compares  the  LS-DYNA  result  at  the 
posterior  sternum  with  the  measured  signal  from  one  of 
the  tests.  The  accelerometer  was  an  Endevco 


2 


7270A-20K,  with  an  amplitude  range  of  20x10^  g’s  and  a 
frequency  range  of  50  kHz  (Endevco  Corp.,  2001).  The 
FE  solution  was  sampled  at  20  ps  intervals  to  be 
consistent  with  the  gauge’s  frequency  response.  A 
general  agreement  between  theory  and  experiment  was 
achieved,  although  oscillations  are  observed  in  the 
computations  that  were  not  present  in  the  experiment. 
Less  capable  accelerometers,  in  terms  of  both  amplitude 
and  frequency  response,  were  employed  at  the  sternum  in 
the  other  two  tests  and  at  the  other  three  locations  in  all 
four  tests.  In  general,  at  the  other  three  locations  the  FE 
signal  exceeded  in  amplitude  the  measured  signal.  See 
Raftenberg  (2003)  for  details.  Conclusions  took  the  form 
of  recommended  improvements  to  the  Kevlar  fabric 
modeling  (pursued  in  section  3)  and  to  the  WSTM  thorax 
model  (pursued  in  section  4),  and  suggested  sources  of 
error  in  the  experiments.  The  latter  include  slippage  of 
the  sutured  gauge  and  inertial  loading  provided  by  the 
gauge. 


time  (ms) 

Fig.  6.  FE  results  and  AFIP  data  for  sternum  acceleration. 


3.  PROGRESS  IN  MODELING 
THE  KEVLAR  VEST 

In  section  2  the  Kevlar  vest  was  modeled  using  isotropic 
linear  elasticity.  Figure  7  provides  a  sketch  of  the  plain- 
weave  construction  of  each  ply  of  the  vest.  Plain-weave 
consists  of  two  initially  mutually  orthogonal  families  of 
yams  called  warp  and  fill,  or  weft.  Note  the  coordinate 
system  defined  in  the  figure.  Material  coordinates  Xi,  X2, 
and  X3  coincide  with  the  warp  yams,  fill  yams,  and  the 
transverse  (through-thickness)  direction,  respectively. 
Based  on  this  constmction,  one  would  expect  material 
orthotropy  rather  than  isotropy  to  characterize  this  vest. 
Figure  8,  reproduced  from  Raftenberg,  et  al.  (in  review), 
presents  stress-strain  data  from  a  single  ply  of  plain- 
woven,  600-denier  Kevlar  KM2  pulled  in  quasi-static, 
uniaxial  tension  along  its  warp  yams.  Figure  9, 


reproduced  from  Raftenberg,  Scheidler,  and  Moy  (2004), 
presents  stress-strain  data  from  the  multi-ply  Kevlar  vest 
subjected  to  quasi-static,  uniaxial,  transverse 
compression.  Figures  8  and  9  make  clear  that  the  vest’s 
stress-strain  response  is  nonlinear  as  well  as  anisotropic. 
The  in-plane  tensile  nonlinearity  in  Fig.  8  can  be 
attributed,  at  least  in  part,  to  progressive  yam  decrimping. 
The  transverse  compressive  nonlinearity  in  Fig.  9  can  be 
at  least  partly  attributed  to  the  progressive  squeezing  out 
of  air  from  the  vest. 

Hence,  ARE  is  engaged  in  the  development  of  a  vest 
constitutive  model  that  is  nonlinear,  orthotropic,  and 
hyperelastic.  The  entire  vest  is  viewed  as  a  homogeneous 
continuum  (Fig.  10).  The  model  is  assigned  the  same 
initial  thickness  Zq  as  measured  on  the  physical  vest 
using  a  micrometer,  an  initial  density  po  given  by 


Fig.  7.  The  plain-weave  constmction  of  a  single  ply,  and 
a  Cartesian  material  coordinate  system. 


Fig.  8.  Second  Piola-Kirchhoff  stress  versus  Green  strain, 
both  along  the  warp-yam  direction,  for  a  ply  of  plain- 
woven,  600-denier  Kevlar  KM2  pulled  in  quasi-static, 
uniaxial  tension  along  the  warp-yam  direction. 
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Fig.  9.  Second  Piola-Kirchhoff  stress  versus  Green  strain, 
both  along  the  transverse  direction,  for  a  multi-ply  vest 
composed  of  plain-woven,  600-denier  Kevlar  KM2  quasi- 
statically,  uniaxially  compressed  along  the  through¬ 
thickness  (transverse)  direction. 


mANSVERSE 


WARP 


Fig.  10.  The  vest  modeled  as  a  homogeneous,  orthotropic 
continuum. 


vest' s  areal  density  .  ^  _ 

Po  = - -J -  >  (1) 

and  we  seek  a  constitutive  relation  of  the  form 
S  =  (p(E).  (2) 

Here,  E  is  the  Green  strain  tensor,  S  is  the  second  Piola- 
Kirchhoff  stress  tensor,  and  (p  is  a  tensor  function.  (See 
Malvern  [1969]  for  definitions  of  E  and  S.) 


Next,  motivated  by  the  plain- weave  construction  in 
Fig.  7,  we  introduce  the  decoupling  assumption,  whereby, 
each  component  of  S  is  a  function  only  of  the  single 
corresponding  component  of  E.  More  specifically,  in 
terms  of  the  coordinate  system  defined  in  Fig.  10, 


^11  - 

(3a) 

^22  ~  ^fi.^22) 

(3b) 

^33  -  (^33  ) 

(3c) 

^12  -  ^12 

(3d) 

S23  =  2G^^23 

(3e) 

A'3i=2G^^31 

(3f) 

This  constitutive  model  introduces  two  constant  shear 
moduli,  G^f  and  G^,  and  three  nonlinear  functions, 
and  Warp  function,  (j)^,  and  transverse  function,  (l)t,  are 
sketched  in  Figs.  11  and  12,  respectively.  Fill  function, 
(l)f,  has  the  same  qualitative  features  as 


Fig.  11.  Warp  function,  (j)^. 


Fig.  12.  Transverse  hmction, 
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(j)^  in  Fig.  11  can  be  described  in  terms  of  four 
subdomains:  compression  (^ii<0),  tension  prior  to 

“rupture”  ( 0  <  i  post-rupture  tensile 

unloading  ( <E,,<  ),  and 

completely  failed  material  (  i  )• 

The  warp  yams  are  assumed  to  buckle  immediately  under 
axial  compression,  so  that  =  0  throughout  the 

compressive  subdomain.  For 

consists  of  a  rational  function  that  introduces  six  material 
constants:  and  f^.  Mathematica  software 

(Wolfram,  1999)  was  used  to  apply  nonlinear  regression 
to  evaluate  these  constants.  The  results,  given  in  Table  1, 
were  found  to  produce  a  good  fit  to  the  uniaxial  tension 

data  throughout  this  subdomain  (Fig.  8).  the 

strain  that  defines  the  boundary  between  the  second  and 
third  subdomains  and  that  corresponds  to  the  maximum 
load  level  reached  in  the  tension  test,  constitutes  the 
seventh  material  constant  associated  with  (j)^.  The  value 

assigned  to  in  Table  1  is  a  direct  measurement. 

For  ^  the  model  assumes  a  subsequent  linear 

ramping  down  of  stress  with  increasing  strain  (Fig.  8). 
This  ramp  introduces  an  eighth  material  constant,  strain 

increment  ,  corresponding  to  completed 

specimen  unloading,  i.e., 

<t>^=0  for  1  >  ).  may  be 

amenable  to  direct  experimental  evaluation,  but  the  value 
of  0.4  given  in  Table  1  was  chosen  to  avoid  large  element 
distortion  in  the  ballistic  simulation. 

Fill  function  ^describes  the  stress-strain  response  in 
a  quasi-static,  uniaxial  test  in  which  a  single-ply  specimen 
is  strained  along  its  fill  yams.  As  was  mentioned,  ^was 
found  to  be  qualitatively  similar  to  the  of  Fig.  11.  ^ 
introduces  the  eight  analogous  material  constants  o/,  Zy,  c/, 

dfi  efi  ffi  ,  and  .  The  first  six  were  again 

evaluated  by  nonlinear  regression  performed  with 
Mathematica  software.  was  determined  directly 

from  the  tension  test.  The  value  of  0.4  used  for 

was  again  chosen  to  avoid  large  element  distortion  in  the 

ballistic  simulation. 

The  transverse  function  in  Fig.  12  is  divided  into 
two  subdomains:  compression  ( - 1/2  <  ^33  <  0 )  and 
tension  ( 0  <  ^33 ).  is  assumed  to  be  zero  in  tension 
because  the  plies  of  the  Kevlar  vest  are  generally  only 
sparsely  stitched  together.  Note  that  -1/2  is  an 
unattainable  lower  bound  on  ^33  that  corresponds  to  the 
state  of  zero  thickness.  In  the  allowable  compressive 


subdomain,  is  governed  by  a  function  introducing  six 
material  constants,  at,  bf,  Ct,  dt,  et,  and  ft.  This  function  has 
a  singularity  at  ^33  =  - 1/2  to  associate  infinite  stress  with 

infinitesimal  thickness.  Constants  at,  bt,  Ct,  dt,  Ct,  and  f 
were  evaluated  by  nonlinear  regression  using 
Mathematica.  The  values  are  added  to  Table  1,  and  the 
resulting  fit  is  seen  in  Fig.  9  to  be  good  over  the  entire 
range  of  compressive  data. 


Table  1.  Material  constants  for  the  multi-ply, 
600-denier,  Kevlar  KM2  vest 


(GPa) 

1.35322 

af  (GPa) 

1.79587 

(GPa) 

-98.8705 

bf  (GPa) 

-204.612 

Cw  (GPa) 

2111.95 

Cf  (GPa) 

14539.8 

d„  (GPa) 

-4898.41 

df  (GPa) 

-33428.5 

-22.0527 

e/ 

-34.3400 

fw 

465.297 

// 

2137.94 

rupture 

^11 

0.132930 

^rupture 

^22 

0.149961 

A  77  unload 
AiSii 

0.4 

0.4 

at  (MPa) 

1.25770 

G„/  (GPa) 

1.0 

bt  (MPa) 

-7.68533 

Gt  (GPa) 

50. 

c,  (MPa) 

-71.1591 

d,  (MPa) 

-135.116 

4.74248 

ft 

6.00453 

Equations  6d  -  6f  relate  each  shear  stress  component 
linearly  to  its  corresponding  shear  strain  component.  This 
introduces  the  constant  elastic  shear  moduli,  G^f  and  Gt. 
The  former  governs  in-plane  shear,  i.e.,  relative  rotation 
between  warp  and  fill  yams.  Currently  we  have  no 
constitutive  data  with  which  to  evaluate  G^f.  We  have 
assigned  it  the  relatively  small  value  of  1.0  GPa  in  order 
to  incorporate  the  “scissoring”  phenomenon,  i.e., 
reasonably  unimpeded  in-plane  relative  rotation  of  the 
two  families  of  yams.  Shear  modulus  Gt  was  assigned  its 
large  value  of  50  GPa  in  order  to  avoid  element  collapse 
in  the  ballistic  simulation. 

Finally,  the  role  of  the  physical  vest’s  multi-ply 
constmction  in  reducing  the  bending  stiffness  must  be 
considered.  According  to  linear  elastic  plate  theory,  a 
plate’s  flexural  rigidity,  which  measures  its  resistance  to 
bending,  increases  as  the  plate’s  thickness  cubed 
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(Timoshenko  and  Woinowsky-Krieger,  1959).  Let  h  be 
the  undeformed  thickness  of  each  ply  of  the  physical  vest, 
and  let  N  be  the  number  of  plies  in  the  vest.  If,  in  the 
physical  vest,  each  ply  behaves  as  an  independent  plate 
with  its  own  neutral  surface,  then  the  net  flexural  rigidity 
of  the  vest  is  KNh^,  where  AT  is  a  constant  determined  by 
material  properties.  On  the  other  hand,  if  one  were  to 
model  the  entire  vest  as  a  single  plate  with  thickness  Nh 
(=Ao),  the  resulting  bending  stiffness  of  the  model  would 
equal  K(Nhf,  in  general  a  much  larger  quantity  than  KNh. 
To  mitigate  excessive  bending  stresses,  we  introduce 
internal  frictionless  contact  surfaces  (“slidelines”)  within 
the  vest.  In  Fig.  13  we  show  four  such  slidelines,  and  the 
vest  is  correspondingly  modeled  with  five  8-node  brick 
elements  through  its  thickness. 


Fig.  14.  Measured  back- face  deflection  of  the  vest 
following  impact  by  an  M882  bullet  at  370  m/s. 


r  (mm)' 


Fig.  13.  Finite  element  models  for  the  vest  and  the  M882 
bullet. 

In  order  to  benchmark  the  FE  fabric  model  prior  to 
placing  it  on  a  thorax,  we  applied  the  model  in  an  LS- 
DYNA  simulation  of  a  ballistic  test  involving  the  vest 
impacted  by  a  M882  at  370  m/s  (Raftenberg,  Scheidler, 
and  Moynihan,  2004).  A  Phantom  V5  camera  was 
employed  to  capture  silhouettes  (2D  projections)  of  the 
vest’s  back  face.  These  stills  were  then  digitized  at  60-ps 
intervals  to  obtain  back-face  deflection  data  (Fig.  14).  In 
the  test  the  bullet  arrested  and  began  to  recoil  at  680  ps 
after  impact.  At  that  time,  the  vest’s  back-face 
displacement  at  the  point  directly  beneath  the  bullet  had 
attained  a  value  of  about  52  mm. 

In  the  LS-DYNA  simulation  the  M882  bullet  was 
approximated  by  a  cylindrical  disc  of  lead  having  an 
initial  diameter,  Dq,  equal  to  the  final  diameter  of  the 
bullet  recovered  after  the  test  (Fig.  15).  The  disc’s  initial 
height,  ho,  was  chosen  to  according  to  the  relation 

=8grams  ,  (4) 


Fig.  15.  The  M882  bullet  was  approximated  by  a  disc. 


By  70  ps  after  impact,  the  time  step  in  the  LS-DYNA 
simulation,  determined  from  the  Courant  stability 
condition,  had  effectively  diminished  to  zero.  The 
solution  up  to  that  time  is  compared  with  experiment  in 
Fig.  16.  Near  agreement  with  experiment  was  achieved  in 
the  rate  of  radial  spread  of  the  displacement  field.  I.e.,  the 
radial  distance  traveled  by  the  Z-displacement  wave 
during  the  40  ps  interval  between  120  and  160  ps  in  the 
experiment  is  about  equal  to  the  computed  radial  distance 
traveled  during  the  30  ps  interval  between  30  and  60  ps. 
However,  there  are  two  notable  points  of  disagreement 
between  theory  and  experiment.  First  is  the  time  required 
for  initial  appreciable  back-face  displacement,  which  was 
between  80  and  120  ps  after  initial  impact  in  the  test  and 
between  10  and  20  ps  in  the  computations.  The  second 
point  of  disagreement  is  in  the  magnitude  of  the  back-face 
displacement.  While  the  peak  magnitude  in  the  test  was 
about  52  mm  (Fig.  14),  at  70  ps  a  magnitude  of  only  11 
mm  was  attained  in  the  computations.  The  bullet  in  the 
computations  was  still  moving  at  70  ps,  but  its  center-of- 
mass  velocity  had  diminished  by  about  78%,  from  its 
initial  value  of  370  to  about  80  m/s  (Fig.  17).  Hence  it  is 
unlikely  that  a  displacement  of  52  mm  would  eventually 
be  attained. 

The  conclusion  is  that  the  FE  fabric  model  is  too  stiff. 
Possible  sources  of  the  excess  stiffness  include  the  neglect 
of  rate  effects,  the  neglect  of  multi-axial  coupling,  and 


where  pie^d  is  the  undeformed  density  of  lead. 
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excess  bending  stiffness  if  an  insufficient  number  of 
internal  slidelines  were  employed. 

2  i  ‘  * . ■  .  .  I  .  . . I  ■  ■  .  >  I 


•14 

Q  5  10  15  a  25  M  35  40 

r  {mrn) 

Fig.  16  The  vest’s  back-face  deflection  following  impact 
by  an  M882  bullet  at  370  m/s;  LS-DYNA  results 
compared  with  experiment. 


Ume(E-03S 


failed  in  that  spurious  penetration  of  internal  contact 
surfaces  occurred. 


Fig.  18.  M80  7. 62-mm  NATO  round. 


To  address  this  challenge,  we  are  currently 
developing  an  FE  thorax  model  that  is  simpler  than 
WSTM.  The  aim  is  to  discard  parts  of  secondary 
structural  importance  and  thereby  create  a  more  robust 
model  having  fewer  contact  surfaces. 

The  revised  thorax  model  is  divided  into  internal 
parts  (Fig.  19)  and  an  external  casing  (Fig.  20).  The 
internal  parts  preserve  WSTM’s  representation  of  the 
circulatory  and  respiratory  systems.  However,  these  parts 
are  now  more  finely  meshed,  with  a  typical  edge  length  of 
5  mm,  reduced  from  10  mm  in  WSTM. 

The  external  casing  consists  of  the  spinal  column,  the 
diaphragm  forming  the  “floor”  of  the  thoracic  cavity,  and 
an  enclosing  shell  part  that  represents  the  ribs  and  the 
muscles  lining  the  ribcage. 


- FE:  10 

-  FE;  Hliia 

-  FE.  30j^S 

FE  Mji* 

FE:  SQjia 

-  FE  SQiia 

-  FE  9e  i|i» 

-  EXPT  120^15 

— - EXFT-.  leoj,* 

- &.fn  KiOns 


Fig.  17.  The  computed  center-of-mass  velocity  of  the 
lead  disc  versus  time. 


4.  A  COMPUTATIONALLY  ROBUST 
THORAX  MODEL 

The  WSTM,  with  its  extensive  anatomical  detail, 
contains  numerous  internal  contact  surfaces  (slidelines) 
between  adjacent  parts.  The  model  functioned  to  produce 
plausible  results  for  the  ballistic  application  of  section  2, 
involving  impact  by  an  M882  at  445  m/s. 

However,  a  second  situation  of  great  interest  is  the 
NIJ  Level  III  design  threat,  namely  an  M80  bullet 
impacting  at  823  m/s.  The  M80  is  a  NATO  round  with  a 
7.62  mm  diameter,  also  composed  of  a  lead  core  and  a 
gilding  jacket  (Fig.  18).  When  WSTM  was  applied  to  this 
more  energetic  impact,  the  LS-DYNA  contact  algorithm 


Ixfi  Lung 


Vena  Cava 
Pulmonary  Anerv' 
Pulmonary  Vein 


Fig.  19.  Internal  parts  of  the  thorax  model. 


7 


‘Thoracic  Wall”  (ribs,  carliiaye.  muscle) 
Abdominal  Muscles 

Vertebrae 


Fig.  20.  The  outer  casing  of  the  thorax  model. 


CONCLUSIONS 

ARL  is  striving  to  develop  an  FE-based  design 
capability  for  thoracic  body  armor.  Towards  this  end, 
progress  has  been  made  in  developing  FE  models  for 
Kevlar  vests  and  for  the  human  thorax. 

Future  plans  regarding  the  fabric  model  include 
numerical  experimentation  to  determine  the  required 
number  of  internal  slidelines.  We  are  also  seeking  to 
expand  the  experimental  database  to  include  multi-axial 
loading  and  high-rate  loading.  Such  new  data  would 
allow  us  to  examine  our  current  assumptions  of 
decoupling  and  rate  independence. 

Future  plans  regarding  the  thorax  model  include  the 
addition  of  internal  fluid  to  represent  the  structural 
contributions  of  cells.  In  addition  we  plan  to  revise  the 
constitutive  modeling  of  various  biomaterials  to 
incorporate  recent  advances  in  the  experimental  database. 
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NIJ  Standard; 

Levels  3A  and  3  Body  Armor 


NIJ  Level  3A  Armor 
o  Design  Threat:  M882  @  425  m/s 

o  Acceptance  Criterion:  44-mm 
Indentation  in  Plastilina  Clay 

o  Usual  Solution:  Kevlar  Vest 

NIJ  Level  3  Armor 
o  Design  Threat:  M80  @  823  m/s 

o  Acceptance  Criterion:  44-mm 
Indentation  into  Plastilina  Clay 

o  Usual  Solution:  Kevlar  Vest  + 
SAPI  Plate 


j/L 
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The  Goal 


Potential  advantages  over  the 
current  laboratory-based  NIJ 
methodology: 

-  Take  into  account  hit 
location. 

-  Less  costly  to  iterate  on 
design. 


Stresses  and 
Displacements 


I - - 

I  Injury 
I  Assessment 


Wayne  State  Thorax  Model  (WSTM99) 


•  Developed  by  Kevin  Wang 
and  King  Yang 

•  Anatomically  detailed  FE 
mesh 

•  Most  materials  represented 
by  isotropic  linear 
elasticity 

•  Typical  element  size  is  10 
mm 

•  Originally  developed  for 
automotive  applications 
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Putting  a  Vest  on  WSTM99 


10-mm  in-plane  element  dimension  compatible 
with  WSTM99 

3.6-mm  initial  thickness  computed  using 
thickness  = _ vest’s  amaldensi^ 


Kevlar's  volumetric  density 

Isotropic  Hooke’s  law  using  Kevlar  KM2  data 
from  Johnson,  Beissel  &  Cunniff  (1999). 


t-f  (-!  1 1  /  T_'. 


LS-DYNA  Simulation: 
M882  vs.  Kevlar  Vest  +  WSTM99 


•  Bullet  modeled  with  1  element  and 
Johnson-Cook  properties  for  lead. 

•  Acceleration  results  were  compared 
with  data  from  AFIP  cadaver  tests. 

•  Some  agreement  at  the  sternum; 
elsewhere  computed  accelerations 
exceeded  measurements. 

•  Conclusions  involved  recommended 
instrumentation  and  refinements  to 
WSTM99  and  to  fabric  modeling. 

•  Documentation: 

Raftenberg,  ARL-TR-2897,  Jan  2003. 


Velocity  Results  from  LS-DYNA  Simulation: 
M882  vs.  Kevlar  Vest  +  WSTM99 


-S.OOOe-OI 

-1.000e.0o1 

Wave  reaches  T4  vertebra  by  traveling  along  4*  rib. 
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Accelerometer  Locations  in  AFIP  Tests 


•  Posterior 
sternum 


T7 

vertebra 


Ligamentum 

arteriosum 


Carina  of 
trachea 
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WSTM99  Contact  Breaks  Down  Under  More 
Energetic  Threats 


1  kg  TNT  @  3  m  standoff 
vs.  WSTM99 


M80  @  823  m/s  vs. 
WSTM99  +  Kevlar  +  SAPI 


Critique  of  WSTM99  As  Applied  to  Ballistics 


•  Coarse  (10  mm)  meshing 

•  No  bending  stiffness  in  ribs 

•  Artificially  large  densities  (low  wave  speeds)  assigned  to  certain 
parts  to  account  for  inertia  of  omitted  anatomy 

•  Vacuum  replaces  fluids  (fat  cells,  blood,  air) 

•  Linear  elasticity  applied  to  most  biomaterials  (neglects  amplitude 
nonlinearity,  rate  effects,  and  damage) 

•  Not  robust  (interface  contact  breaks  down) 


6 


w 


Suggested  Simplifications/Improvements 
to  the  WSTM 


Two  criteria  for  what  to  keep: 
o  If  you  don’t  keep  it,  you  won’t  get  output  for  it. 

-  Heart  T 

-  Lungs  >  sites  of  “not  immediately  survivable  ’’  injuries 

-  Liver  J 

o  Does  its  presence  affect  the  solution  for  the  parts  of  interest? 

-  Associated  Plumbing  -  anchor  the  Heart  &  Lungs 

-  Fat,  Blood  &  Air  -  transmit  load  to  Heart  &  Lungs 

-  Ribcage  and  the  Muscles  that  line  it  -  main  load-bearing  elements 

-  Diaphragm 

Spine  thoracic  cavity 


Outer  Casing 


X 

Left  view 


Front  view 


New  material 
represents 
combination  of 
ribs  and  lining 
muscles  -  to  be 
tuned  to  UVa 
cadaveric  data 


Refined  FE  Model  Development 
for  a  Fabric  Vest 


ARL  co-workers: 

Theoretical/computational:  Michael  Scheidler 

Experimental:  Thomas  Moynihan,  Thomas  Mulkem,  Paul  Moy 

•  We  seek  a  constitutive  model  for  the  vest  that  is  useable  in  large- 
scale  ballistics  calculations  involving  a  bullet,  a  vest,  and  a  human 
torso. 

•  In  the  interest  of  computational  efficiency,  we  would  prefer  that  the 
fabric  vest  be  represented  as  a  homogeneous  continuum. 


vest’s 

thickness 


A 


vest’s  areal  density 
vest’s  initial  thiekness 

S  =  (p(E) 
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-0.5<E„<0 


2S  Plies  erf  Plain-Woven,  SOO-Denier,  Kevlar  KM2  {Style  706) 


2B  Plies  of  Plain-Woven,  600-Denier,  Kevlar  KM2  (Sty  le  706) 


Material  Constants 


«w(GPa) 

1.35322 

^w(GPa) 

-98.8705 

c„(GPa) 

2727.95 

^w(GPa) 

-4898.41 

-22.0527 

f 

Jw 

465.297 

r  rupture 

^ww 

0.132930 

a^(GPa) 

1.79587 

Zy(GPa) 

-204.612 

^(GPa) 

14539.8 

^.(GPa) 

-33428.5 

G- 

-34.3400 

// 

2137.94 

r  rupture 

0.149961 

(MPa) 

1.25770 

b,  (MPa) 

-7.68533 

c,  (MPa) 

-71.1591 

d,  (MPa) 

-135.116 

G 

4.74248 

f, 

6.00453 

G,,^(GPa) 

1.0 

G,  (GPa) 

50. 
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Applied  individually  to  each  non-eroded  fabric  element  at  each 
time  step: 

Tf  F  >  C’WPture  A  p  unload 

—  ^ww--^ww 

thereafter  the  element  can  no  longer  support 
-If 

thereafter  the  element  can  no  longer  support 

-  If  both  conditions  have  been  met,  the  element  is  eroded 
from  the  mesh. 


A  Meshing  Consideration; 
Mitigation  of  Bending  Stresses 


h  =  ply  thickness 

Bending  stiffness  of  vest 
with  A  plies  az  Nh^ 

Bending  stiffness  of  a 
single-ply  vest  with  the 
thickness  Nh  oc  {Nhf 


4  frictionless  slidelines 
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Multi-axial  Effects? 


Quasi-static  multi-axial  testing  is 
needed  to  check  the  decoupling 
assumption. 

-  Biaxial  in-plane  tension 
tests  have  been  performed 
on  other  fabrics. 

-  Tri-axial  tests 
(superimposed  transverse 
compression)  have  not  yet 
been  done. 


From  Boisse  et  al,  NSF  Workshop  on  Composite  Sheet 
Forming,  Sept.  5-7,2001. 
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Backface  Deflection; 

M882  @  380  m/s  vs.  Fabric  Vest 


•  Visual  Solutions  V-5  Camera 


Backface  Deflection  Data  From  the  Movie 
(M882  @  370  m/s  vs.  Fabric  Vest) 


-  120  tis 

- 160  tis 

-  200  |is 

-  240  trs 

-  280  |is 

-  320  tis 

-  360^15 

-  400  |is 

-  440  jis 

-  480  tis 

-  520  tis 

-  560  trs 

-  600  tiS 

-  640  iis 

-  680  IIS 


-5 
-10 
-15 
-20 
-25 
N  -30 
-35 
-40 
-45 
-50 
-55 
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FE  Shape  Function  Suppresses  Mushrooming 


Vx  y)  =  ax+  b^x  +  c^y  +  d^xy 


V;,  (x,  y)  =  ay+  byX  +  j  +  dyxy 


•  As  the  element  stretches  over  a  large  region,  its  low-order 
interpolation  functions  introduce  spurious  stiffening. 

•  Three  standard  approaches  to  rectify  this 

-  Automatic  remeshing 

-  Particle  methods  such  as  SPH 

-  Lagrangian-Eulerian  hybrid  such  as  ALE 
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Backface  Deflection:  FE  and  Experimental  Results 
(Disc  @  370  m/s  vs.  Vest) 


FE.  lOpis 

PIT*  'XI 

rc.  ^jis 

f-f-.  0%r\ 

rt.  oU  )1S 

rc.  4U 

CC*  lie 

rc.  ou  jis 

rt.  tSU  jiS 

FE.  70 

FE.  80  ixs 

f\r\  . .  ^ 

rt.  w  JiS 

rt.  yo  JiS 

t  Ar  1  .  1  |1S 

- EXPT;  160  ns 

t  An  1 .  zUU  |iS 

tAn  1  .  ^4U  |iS 

- EXPT:  280  fis 
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Fabric  Modeling: 

Error  Analysis/Planned  Improvements 


•  Features  of  the  Material  Model 

-  No  rate  effects  1  Awaits 

-  No  multi-axial  coupling  effects  J  new  data 

•  Features  of  the  FE  Application 

-  Pre-mushrooming  the  bullet  (plastic  work,  the  bullet’s 
early-time  shape)  -  SPH  or  ALE  should  work 

-  Spurious  bending  stiffness  if  not  enough  slidelines  - 
numerical  experiments  will  study  this 

-  Lateral  dimensions  of  target  must  be  increased 
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